A Cosmological Upper Bound on Superpartner Masses 
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If some superpartners were in thermal equilibrium in the early universe, and if the lightest su- 
perpartner is a cosmologically stable gravitino, then there is a powerful upper bound on the scale 
of the superpartner masses. Typically the bound is below tens of TeV, often much lower, and has 
similar parametrics to the WIMP miracle. 



INTRODUCTION. A natural weak scale, precision 
gauge coupling unification and dark matter provide pow- 
erful arguments for weak-scale supersymmetry. However, 
to date, direct evidence for supersymmetry is still miss- 
ing and thus whether or not low-scale supersymmetry is 
realized in nature remains unknown. In fact, the recent 
discovery of a 125 GeV Higgs boson [TJ [5] implies a fine- 
tuning in the MSSM worse than 1% [5], and searches for 
supersymmetry at the LHC are placing limits on colored 
superpartners in the region of 1 TeV [HE]. Therefore, the 
superpartner mass scale, m, may be decoupled from the 
weak scale and could, in principle, be anywhere between 
the present experimental bound near 1 TeV up to the 
Planck scale. With the naturalness reasoning aside, the 
question arises: Are there arguments for superpartners at 
the TeV scale that are unrelated to the stabilization of the 
weak scale? 

The argument for TeV superpartners from gauge cou- 
pling unification alone is weak, as logarithmic running 
implies that the precision changes only mildly as m in- 
creases well above 1 TeV. On the other hand, there is a 
powerful and well-known argument for TeV-scale super- 
partners from the cosmological abundance of the lightest 
supersymmetric particle (LSP) [6 . This results from LSP 
freeze-out and follows from three assumptions: 

(i) . The LSP is cosmologically stable. 

(ii) . The reheat temperature of visible particles after in- 

flation, Tr, was sufficiently high, Tr > fh. 

(iii) . There is no substantial late-time dilution of the 

LSP abundance. 

The second assumption implies that the standard model 
superpartners were in thermal equilibrium. If we further 
assume 

(iv-A). The LSP reached thermal equilibrium, 

then the thermal freeze-out relic abundance leads to the 
overclosure bound on the LSP mass, 

m\ SP < 2.0^ a 2 cS T eq M P1 ~ (2.3 TeV ^g)" , (1) 

where T cq ~ 1.5 eV is the temperature of matter- 
radiation equality, M P1 ~ 2.4 x 10 18 GeV is the re- 



duced Planck mass, and Xt is fh divided by the freeze- 
out temperature. The coupling strength a c s, appearing 
in the thermally averaged LSP annihilation cross section 
at freeze-out, is defined by (cry) = 47ra 2 ff /m 2 . It is 0.03 
(0.01) for wino (Higgsino) LSP, but can have a much 
larger variation. 

The TeV scale from freeze-out thus results paramet- 
rically as the geometric mean of T eq and Mpi and is 
independent of the weak scale. The equality holds for 
LSP dark matter. Eq. is a very important result: in 
models where superpartner masses are characterized by 
a single scale, fh is likely in the 1-10 TeV window, and in 
Split Supersymmetry [7] the fermonic superpartners lie 
in the TeV region. 

From the above list of assumptions it is clear how to 
evade the bound in Eq. ([T]), allowing fh many orders of 
magnitude above the TeV scale. In particular, violat- 
ing assumption (i) through, e.g., R parity breaking, may 
void the bound entirely. In that case, DM could arise, 
for instance, from a hidden sector or from axions. Violat- 
ing assumption (ii), having m well above Tr, allows the 
superpartners to have no cosmological role, hence evad- 
ing the bound. Finally, assumption (iii) may not hold if 
additional late-decaying states reheated the universe. 

In this paper we study the intriguing possibility of vi- 
olating assumption (iv-A). Indeed, there are numerous 
scenarios where DM is only very weakly coupled so that 
its abundance does not follow from thermal freeze-out, 
invalidating Eq. ([IJ. 

The most common scenario of this kind has (iv-A) re- 
placed by 

(iv-B). The gravitino is the LSP (and the Lightest 
Observable-sector SuperPartner (LOSP) decays 
predominantly to gravitinos). 

The gravitino, present in all supersymmetric theories, has 
interactions that are highly constrained and very weak. 
The gravitino has a cosmological abundance determined 
by thermal scattering, freeze-in, and freeze-out and de- 
cay, and reaches thermal equilibrium, in accordance with 
(iv-A), only when it is very light. The gravitino abun- 
dance has been studied in detail for the case of weak- 
scale superpartners, for example leading to bounds on 
Tr as a function of the gravitino mass [S]. In this let- 



2 



Bound on Superparticle Mass Scale in Single Scale SUSY 



io 5 



10- 



> 

O 



10- 



An 



gsusy" = 10 



overdosed / 

/ 

/ 



<*eff = 0.03 (wino) 



Th 



/ 



< — therraalized 
G / 
/ 



Log,, 



/' ,-' /'TJW 



10 



Ssusy ~ 3 / £i 



£susy — 1 



/ 



/ 



FI, 



/ / / ,< 
/ / / / 



, 0/ .1 2 3 4 .5 6 .7 / 




10 -13 10 -12j -ll 10 -10 1Q -9 j Q-8 1Q -7 j q-6 10 -5 jq-4 jq-3 1q -2 jq-1 1 jqI 1Q 2 10 3 10 4 10 5 

m 3/2 [GeV] 



FIG. 1: The cosmologically allowed region in the (1713/2, fh) plane, for a single scale SUSY with an LSP adhering to assumptions 
(i)-(iv-B) discussed in the text. The gravitinos are (are not) thermalized to the left (right) of the orange dashed line (assuming 
Tr = m). Even when Tr is as low as fh, gravitinos provide too much dark matter in the red region, which has borders labelled 
by the relevant process Th, FI or FO. As Tr is increased the overdosed region becomes larger, as illustrated by the dashed 
blue lines, because UV scattering at Tr produces more gravitinos than freeze-in. At the edge of the red region (suitably 
enlarged for Tr > fh) gravitinos provide the observed dark matter. In the region to the right of the slanted black dashed line 
the gravitino is not the LSP; this is the conventional WIMP LSP freeze-out region, with a limit of 2.3 TeV for a wino LSP. 
The green region is excluded by the effects of late decays of LOSPs to gravitinos during big bang nucleosynthesis (BBN) [18] ; 
light green shading corresponds to a neutral LOSP with 100% hadronic BR, and dark green shading to a neutral LOSP 
with 1% hadronic BR and 99% electromagnetic BR. The BBN limits when the gravitino is not the LSP are model dependent 
and are not shown [T5] . The purple region next to the "Th" contour is excluded as the gravitino component of dark matter is 
too warm [20j . The gray shading (and corresponding gray dashed and dotted lines) shows the regions with gl usy > 10 (3, 1), 
which are excluded as described in the text. 



ter, however, we take a different approach and derive the 
cosmological bound on the superpartner mass scale for a 
gravitino LSP. We find this bound to be strong, so that 
under the quite mild assumptions of (i), (ii), (iii) and 
(iv-A) or (iv-B), supersymmetry, if it exists, must be in 
the (multi-) TeV domain. We also derive bounds for the 
split spectrum case and scenarios where the LOSP does 
not predominantly decay to the gravitino. 
SINGLE SCALE SUSY. In this section we take all su- 
perpartners of the observable sector to be characterized 
by a single mass scale, fh, and leave the case of a non- 
degenerate spectrum to the next sections. Our aim is to 
derive a general bound on the scale fh from overproduc- 
tion of gravitinos. We ignore other possible components 
to DM since they would only strengthen the bound. A 
key superpartner is the LOSP, since it undergoes freeze- 
out. We allow a very wide variation in the (TO3/2, in, Tr) 
space. 

The upper bound on fh follows from the three as- 
sumptions (i), (ii) and (iii). Assumption (ii) implies 
that the observable sector produces gravitinos from three 
sources: gaugino scattering at Tr [SHU]: ¥3/21 gravitino 
"freeze-in" from decays of visible sector superpartners 



at T ~ m [13], Yfl, and LOSP freeze-out and de- 
cay |14) . Yy® ■ F° r sufficiently small m 3 / 2 , the grav- 
itinos are in thermal equilibrium when T = fh; in this 
case Y^, v + Yy^ are replaced by a thermal abundance, 
and ma Y be neglected. Below, in accordance with 

assumption (iv-B), we assume the LOSP branching ra- 
tio to the gravitino is 0(1). In the final section we 
discuss how our bound is weakened when this assump- 
tion is relaxed. Gravitinos may also be produced from 
other sectors or they may arise from an initial condi- 
tion [121 [H]. However, these additional sources of grav- 
itinos only strengthen our bound, and to be conservative 
we ignore them. 

If gravitinos do not thermalize, the condition that they 
not yield too large a DM abundance is 

Trtti 2 fh 3 mm 3 / 2 
Cuv 1- Cfi \~ Cpo 2 — — o,Mp\T eq , 

m 3 / 2 TO3/2 Qfeff 

(2) 

where a — 0.27 and a e g is now the coupling relevant for 
LOSP annihilation. The three terms labelled UV, FI and 
FO correspond to scattering at Tr, freeze-in and freeze- 
out and decay and occur with rate constants Cuv = 
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FIG. 2: Left: The bound on rh in the single-scale SUSY case, for a e ff = 0.03, 10 -2 and 10~ 3 in blue, green and purple 
respectively, assuming Tr = rh. As a e ff decreases freeze-out yields a larger abundance, so the FO boundary and the BBN 
constraints (shown shaded in the corresponding colors) both become more stringent. As Tr is raised, the bounds become more 
stringent as indicated by the blue dashed lines of Fig. [I] Right: The bound on rh when the contribution to the gravitino 
abundance from freeze-out and decay is negligible. This may be the case in several scenarios, as discussed in the last section. 
The dashed blue lines demonstrate the strengthening of the bound as Tr is increased. We do not analyze the region with 
rh < m 3 /2 as the results are model-dependent. 



73 2^3/2 ^ 2.4 x 10 , L F i - ^y^-^m^r - 

3.8 x 1CT 4 and C FO = 8 ^>2 ~ O- 13 ^)- Here 
73 ~ 0.36 is related to the thermal corrections of the 
scattering process [TT], = 228.75, and npi counts the 
number of fermions and complex scalars participating in 
the freeze-in with mass m; with degenerate MSSM spar- 
ticles, n FI = 36 + 9 + 12 + 4 = 61. The equality in Eq. Q 
corresponds to the case that these processes yield the ob- 
served DM abundance. If gravitinos do thermalize, the 
overabundance constraint becomes [17] 

Cth %/2 < aT cq , (3) 

with C Th =Yy = 45£(3)/7r 4 3>ts w 2.4 x 10~ 3 . Here g* s ~ 
<?* = 228.75. The resulting bound on rh as a function 
of m 3 / 2 is shown in Fig. fTlfor a c fr = 0.03, relevant for a 
(perturbative) wino LOSF. We do not include the non- 
perturbative Sommerfeld effect [21j . which results in an 
O(l) shift in a c ff. 

When gravitinos are not thermalized, the key point is 
the differing dependences of the three terms in Eq. ([2| on 
rh and 7713/2- While all three terms have a positive power 
of rh, the UV and FI terms are proportional to 1/7713/2 
while the FO term is proportional to 7773/2, leading to 
contours in Fig. [T] with slopes of opposite signs. Hence 
there is an upper bound, 

m 2 < ^ 2 a cS Mpi T eq , (4) 
vCfo^d 



where Co — Cuv(TR/rh) + Cfi- At the bound m 3 / 2 = 
\/Cd/Cfo a cS 777,. For Tr 3> 771 the bound becomes rh < 
27TeV[(T i? /?77)/10]- 1 / 4 for a cff = 0.03 which weakens 
to 777, < 38 TeV for Tr = rh. Decreasing a e g makes 
the FO term larger, as shown in the left panel of Fig. [2] 
for Tp; = rh. The parametrics of Eq. |4]) is similar, but 
not identical, to that in the so-called "WIMP Miracle" , 
Eq. 0. 

A second allowed region occurs at very low 777,3/2 in 
Fig. [I] where the gravitinos are thermalized for any 
Tr > 771. Here the bound on 777 arises from theory rather 
than cosmology: rh < (g S usy/^) 2 vF, where g SU sy is the 
strength of the coupling between obervable and super- 
symmetry breaking sectors, and F = ^/im^^M-py is the 
supersymmetry breaking scale. The bound results when 
the messenger scale takes its minimal value of V~F, and 
is shown in Fig. [I] for g 2 usy = 1,3 and 10. We note 
that it may be possible to construct realistic models of 
composite quarks and leptons having non-perturbative 
couplings, g susy ~ 4tt [22]. 

NON-DEGENERATE SPECTRUM. The com- 
pletely degenerate spectrum discussed above is special 
because non-degeneracies typically arise from renormal- 
ization group effects or the dynamics of the mediation 
of supersymmetry breaking. How do non-degeneracies 
affect the above bounds? 

Non-degeneracies induce independent changes in the 
three gravitino production mechanisms. The freeze-in 
process is dominated by the heaviest superpartners, ?ti + , 
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FIG. 3: Left: Bounds in the (m 3 y 2 ,mnc) plane for colored (non-colored) superpartners with mass m c (m ne ). The importance 
of freeze-in as m c /m nc is raised from 1 to 10 is seen by comparing the orange and blue lines. The solid and dashed lines show 
the effect of increasing Tr by 100. Center: Similar to the left panel, the changes to the bound of Fig.[T]is shown for the case of 
split-SUSY, where the scalar superpartner masses, ih s , are raised above the fermionic superpartner masses, rhf. Right: The 
overclosure bound in the (rhs, rhf) plane is shown for the split-SUSY case, where the gravitino mass has been chosen at each 
point to maximize the allowed region. For split-SUSY Tr = rh s . In all panels the green shading is as in Fig. [I] 



and is suppressed compared to the degenerate case by 
n Fi/ n Fi, where ripr is the number of these heavy super- 
partners. The scattering process, dominated by gluino 
scattering, is proportional to the square of the gluino 
mass, Af|. Finally, the freeze-out abundance is propor- 
tional to the LOSP mass, to_, with (av) = 47ra 2 ;ff /m 2 _, 
so that Eq. Q becomes 

T R Mj C FI n+j rh\ m_m 3/2 

Ul/V 1 r^FO 9 S d-Mpi-teq- 

m 3/2 n FI m 3/2 o£ ff 

(5) 

While pure FO of Eq. (JlJ bounds ttilsp, with a gravitino 
LSP the bound depends on the ttilosp, ^3, and the mass 
dominating FI. 

As a simple example, on the left of Fig. [3] we show the 
bound that results by taking all colored states at rh c — 
m + and all non-colored states at rh nc = rh- , assuming all 
superpartners are reheated. As can be seen, the bound 
on rh nc becomes much more stringent as fh c is raised, 
being reduced to 7 TeV for rh c /rh nc — 10. Much of the 
allowed regions in Figs. [T] [2]-Left and [3]Teft are within 
the LHC reach. 

SPLIT SUSY. In the split-SUSY scenario 7 , where the 
scalar superpartner mass, rh s , becomes much larger than 
the fermionic superpartner mass, rhf, a bound on rhf, 
with a gravitino LSP, was discussed in [23]. The freeze- 
in process dominates over the scattering process as long 
as T R > rh s [231 121] • Using Eq. with rh s = rh + 
and rhf = m_, yields the bound on rhf shown in the 
center panel of Fig. [3] for various values of rh s /rhf. To 
compute the bound, the split-SUSY 1-loop RGEs were 
used [25J [25]. The bound on fh s is in the region of 100 
TeV, as shown in the right panel of Fig. [3j and hence 
arbitrary flavor and CP violation in the squark mass ma- 



trix requires Tr < fh s . Finally, we note that if Tr is 
indeed below rh s a bound on rhf may still be obtained, 
and is similar to that shown in Fig. [T] up to 0(1) cor- 
rections stemming from the absence of some diagrams in 
the finite-temperature thermal production of the graviti- 
nos [S]. 

The non-degeneracies explored in the left and center 
panels of Fig. [3] lead to similar bounds, and forbid large 
splittings between the light and heavy states (assuming 
that both are reheated). Indeed, as the splittings in- 
crease, the BBN bounds rapidly become very constrain- 
ing. 

RELAXING ASSUMPTION (iv-B). We now con- 
sider how the bound on superparticle masses is relaxed 
in theories that violate assumption (iv-B). 

LOSP freeze-out and decay may not produce a sig- 
nificant yield of LSP gravitinos, depleting Y^S ■ This 
occurs, for example, if the LOSP dominantly decays 
through i?-parity violating (RPV) operators, which can 
still be consistent with gravitino DM for sufficiently small 
RPV [27J [25]. Alternatively, the LOSP may dominantly 
decay to a light hidden sector, which, if thermalized, may 
not produce significant gravitinos due to its lighter mass 
scale. A third possibility is that the LOSP is colored, in 
which case a late annihilation stage, after the QCD phase 
transition, can dilute the abundance of i?-hadrons [2^1150] 
before the LOSP decays to gravitinos. In these cases, a 
bound on rh results from dropping the FO term and is 
shown on the right of Fig. [2| The maximal rh occurs at 
TO 3/2 — when Eq. ([2]) gives 

m 2 < 7^ T eq M pl < (10 3 TeV) 2 . (6) 

The numerical value above was obtained for Tr = rh. 
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FIG. 4: A bound on rh nc for a sneutrino LOSP with gravity- 
mediated supersymmetry breaking and m c /rh nc — 6 and var- 
ious Tr. The dashed blue lines show the bound for different 
values of Tr normalized to rh c . The purple line at the bottom- 
right corner shows the bound for Tr = 10 9 GeV, correspond- 
ing to the rough reheat temperature required for successful 
thermal Leptogenesis. 



inantly by rare sneutrino decays with a radiated Z or 
W, so the excluded green region is quite small [HI 151] . 
allowing various possibilities. One has a light, e.g. 200 
GeV, sneutrino, with M* near Mpi and a high T fl ~ 10 9 
GeV, compatible with Leptogenesis [35]. In this case the 
colored superpartners may well be in reach of the LHC. 
Another possibility has M* further from Mpi and a much 
lower Tr so that the sneutrino mass can be near its upper 
bound of 5 TeV. 
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For larger reheat temperatures the bound is stronger. 

A more drastic possibility is to consider an LSP that 
violates both assumptions (iv-A,B) entirely, i.e. a state 
that is not the gravitino and yet interacts with the ob- 
servable sector so weakly that it remains out of equilib- 
rium. An example is a light, weakly coupled singlino. In 
this case, the bound can be completely removed as the 
singlino couplings can be chosen to be arbitrarily small 
(removing scattering and freeze-in production) simulta- 
neously with a vanishing mass (thereby removing freeze- 
out and decay), allowing arbitrarily heavy superpartners. 
The key characteristic about the gravitino that leads to 
our bound is that its mass is inversely related to its cou- 
pling to observable states, so that the mass and coupling 
cannot simultaneously be taken too small. 
GRAVITY MEDIATION. When mediation of super- 
symmetry breaking occurs at a very high fundamental 
scale, M*, of order the scale of gauge coupling unifica- 
tion or higher, then m 3 / 2 /m ~ M*/Mp; ~ 10~ 3 — 1. 
Thus "gravity mediation" typically has a gravitino LSP 
and selects a small region of Fig. [T] that is within a few 
orders of magnitude of the m 3 / 2 = rh dashed line. Part 
of this region, with near Mp; , is typically highly con- 
strained by BBN, but smaller values of M* are of interest 
and include the largest values of rh. 

The details of this gravity-mediated region are highly 
dependent on the LOSP, the superpartner spectrum and 
Tr. In Fig. [4] we show a particular example: a sneu- 
trino LOSP with rh c / rh nc = 6. BBN is affected dom- 
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